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We study the electronic transport through a noncollinear single-electron spin- valve transistor. 
It consists of a small metallic island weakly coupled to two ferromagnetic leads with noncollinear 
magnetization directions. The electric current is influenced by Coulomb charging and by spin 
accumulation. Furthermore, the interplay of Coulomb interaction and tunnel coupling to spin- 
polarized leads yields a many-body exchange field, in which the accumulated island spin precesses. 
We analyze the effects of this exchange field in both the linear and nonlinear transport regime. In 
particular, we find that the exchange field can give rise to a high sensitivity of the island's spin 
orientation on the gate voltage. 
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I. INTRODUCTION 

The intense investigations in the field of spintron- 
ics have been partially driven by the prospect of de- 
veloping new devices for current and future informa- 
tion technology*! — Continuing the trend of miniatur- 
ization down to the nanometer scale leads to devices 
in which Coulomb- interaction effects become important. 
Prominent examples are Coulomb blockade and Coulomb 
oscillations. 5 Hence, a fundamental understanding of the 
interplay between spin degrees of freedom and Coulomb 
interaction in nanoscale devices is required. A convenient 
model system to study this interplay is a single-electron 
spin-valve transistor, see Fig. [TJ It consists of a metal- 
lic island that is tunnel coupled to ferromagnetic source 
and drain leads. The current flowing through the island 
is controlled by both an applied gate voltage and the rela- 
tive orientation of the magnetization directions of source 
and drain lead. Experimental measurements in vari- 
ous versions of the single-electron spin- valve transistors 
with collinearly (parallel or antiparallel) polarized leads 
addressed current- volt age characteristics, tunnel magne- 
toresistance (TMR), spin accumulation, and magneto- 
Coulomb effects. This includes devices in which the 
central island is ferromagnetic (FFF)^r— non-magnetic 
(FNF)ji£~— or even superconducting (FSF)j21~— Many 
theoretical works also focused on collinear configurations 
of the leads' magnetization^"— in order to discuss, e.g., 
an oscillating TMR with variation of bias voltage, an en- 
hanced TMR in the Coulomb blockade regime, or a finite 
spin accumulation on the central island. But also non- 
collinear single-electron spin valve transistors have been 
discussed in the literature *22r— 

The single-electron spin-valve transistor is conceptu- 
ally similar to a quantum-dot spin valve. In the lat- 
ter, the central island hosting a continuum of single- 
particle energy levels is replaced by a quantum dot with 
a discrete level spectrum. This can be realized by ei- 
ther shrinking the central island in size or by using 
semiconductors or carbon nanotubes instead of metals. 
Quantum-dot spin valves have been studied extensively 
both theoretically 4 ^"— and experimentally^ 5 -"— One in- 



triguing prediction 4 ^ — for the quantum-dot spin valve 
was the existence of an interaction-induced exchange field 
that acts on the spins of the quantum-dot electrons as 
a consequence of the tunnel coupling to spin-polarized 
leads. This exchange field, which is tunable by gate and 
bias voltage, leads to a precession of an accumulated dot 
spi n 46 i 47 or a splitting in the Kondo resonance^ The 
latter has been experimentally confirmed recently j 67 i 70 i 71 

The origin of the exchange field is a level renormal- 
ization that is spin dependent as a consequence of spin- 
dependent tunnel couplings. This idea has later been 
transferred to describe gate-dependent tunneling-induced 
level shifts in carbon nanotubes with orbital-dependent 
tunnel couplings to normal leads^ and to molecular sys- 
tems with different tunnel couplings of degenerate molec- 
ular states to two normal leads. 75 It is quite natural to ex- 
pect the existence of a similar exchange field for a single- 
electron spin- valve transistor. The implications of such 
an exchange field on the linear conductance has already 
been discussed in the Coulomb-blockade regime j 42 i 43 In 
the present paper, we consider the same system and 
derive kinetic equations for the island charge and spin 
within a diagrammatic real-time transport formalism 
up to the lowest order in the tunnel-coupling strength 
(sequential-tunneling limit). The resulting kinetic equa- 
tions, which agree with Refs. 42.43 in the linear-response 
regime, are the basis for the discussion of the spin accu- 
mulation and its implication on transport in linear and 
nonlinear response. The used formalism systematically 
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FIG. 1: (Color online) Noncollinear single-electron spin- valve 
transistor: a metallic island is tunnel coupled to ferromagnetic 
source and drain leads, whose polarizations directions enclose 
an angle (f>. 
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includes all contributions to a given order in the pertur- 
bation expansion. In addition to the spin-dependent tun- 
neling rates, the exchange field automatically emerges in 
the kinetic equations as a result of the considered Hamil- 
tonian. 



II. MODEL 

The single-electron spin-valve transistor shown in 
Fig. [D is modeled by the total Hamiltonian H = Hi + 
H C + H L + H R + H T . The first part, 



(1) 



describes the metallic island whose energy spectrum e\ 
is characterized by a small level spacing As. The index 
/ labels the energy levels of the island, a G {t, 1} the 
spin, and v = 1, . . . , N c is the transverse channel index. 
The annihilation (creation) operator of island electrons in 

the state lav is denoted by cfj y . We assume that the en- 
ergy levels are independent of spin and transverse channel 
number. We will consider the limit of temperature and 
bias voltage larger than the level spacing, /c#T, eV ^> As, 
for which the energy spectrum can be viewed as contin- 
uous. 

The Coulomb interaction of the island electrons is ac- 
counted for by the charging-energy term 



H c = E c {N-N wA f 



(2) 



where N is the number of electrons on the island. The 
charging energy scale Eg = e 2 /(2Ct) is determined by 
the total capacitance Ce , which is the sum of the capaci- 
tances of the two tunnel junctions, Cx, Cr, and the gate, 
Cg- For equal capacitances of the two tunnel junctions 
and a symmetrically applied transport voltage, the exter- 
nal charge eiV ex t = CqVg depends on the gate voltage 
Vg only. For later convenience, we define An as the dif- 
ference of charging energies of N + 1 and N electrons, 
i.e., A N = E c [2(N-N ext ) + i\. 

Each of the leads is described as a reservoir of nonin- 
ter acting fermions: 
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(3) 



with indices for lead r G {L, R}, momentum k and a^ sv 
is the annihilation (creation) operator of lead r. The 
index s = +(— ) denotes the majority (minority) spin 
states (along the magnetization direction n r of the re- 
spective lead) with the density of states p r s , which we 
assume to be energy independent. The lead's degree of 
spin polarization (at the Fermi energy) is characterized 
by p r = (p r + — p r -) / (p r + + p r J) • The angle enclosed by 
and i\R is denoted by <fi. 

The tunneling Hamiltonian Ht = £ r Ht,v , with 
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FIG. 2: (Color online) Scheme of relation between defined an- 
gles a, /3, the polarization directions of the two leads n^n^, 
and the accumulated spin on the island S. 



describes tunneling between island and leads. Both the 
spin and the transverse channel index v are conserved 
during tunneling. The latter is obvious from the fact 
that the tunneling Hamiltonian is diagonal in v. In 
the following, we assume the tunneling matrix elements 
Vksav = Vsa to be independent of momentum h and 
transverse channel index v. Spin conservation is ac- 
counted for by expressing as a product of the (spin- 
independent) tunnel amplitudes t r and the matrix ele- 
ments of an SU(2) rotation that connects the (in general) 
different spin quantization axes for the two leads and the 
island. For our derivation of the kinetic equations it is 
convenient to choose the spin quantization axis of the 
island along the direction of the accumulated island 
spin S. Its orientation relative to the lead's magneti- 
zation directions can be parametrized by two angles, as 
shown in Fig. [2j the angle a enclosed by the n^- n^-axis 
and the projection of S onto the (n^, n^-plane, and the 
angle f3 between S and the (n^, n^-plane itself. Then, 
the tunnel matrix elements Vg a for the left lead become 



= t L 
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V2 



(5) 



±e i<t>/2 sin (f - f ) + ie ia cos 
V2 



£ _ ZL 



. (6) 



while the elements of the right lead are described by the 
same expressions but with the replacements L — >• R and 
(/>—)>—(/>. The tunneling rate for electrons from lead 
r with spin s into the island spin state a is quantified 
by T r sa /H = 27Tp r s \V s r a \ 2 /H. In addition, we define T r a = 

E s K*, r r = r ; A as well as r = £ r r r . 



III. METHOD 

The dynamics of a quantum-mechanical system is de- 
termined by the time evolution of its total density matrix. 
Since we do not want to investigate the behavior of the 
lead degrees of freedom, we integrate them out and ob- 
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tain a reduced density matrix p re d, which only contains 
the degrees of freedom of the island. An island state 
is characterized by the ket vector \x) = \{ni au }), where 
niav = 0, 1 counts the occupation of the corresponding 
island level. In this notation, the elements of the reduced 
density matrix are defined as P^ = (xi| /W \X2)- Diag- 
onal elements P* = P x describe the probabilities to find 
the island in the corresponding state For later conve- 
nience, we introduce the notation \x) = Kxz^; where 
Xiv describes the occupation of all island levels except for 
the orbital / and channel index ^, while ip G {0,t,i, d} 
explicitly shows the occupation of level / and channel v. 

The time evolution of p re d is determined by the follow- 
ing kinetic equations: 

±P£(t)= - l -(E Xl -E X2 )P£(t) 

+ f dt'j2w*;£(t,t')p£(t').(7) 

>' — OG / / 

XiX 2 

Here, the definition of the energies of the reduced system, 
E x = ( X \ (Hi+Hc) \x), is used. The kernels (t, t' ) 

describe transitions between the matrix elements P Xl (t f ) 

x 2 v ' 

and P x ^(t) of the initial and final density matrix, re- 
spectively. They can be obtained using a diagrammatic 
real-time technique formulated on the Keldysh contour. 
It allows for a systematic perturbative expansion in the 
coupling strength. In this paper, we truncate the expan- 
sion at the lowest order T to describe the weak-coupling 
limit (sequential tunneling). In general, Eq. ([7]) describes 
non-Markovian behavior, i. e., the time evolution of ma- 
trix element Pg (t) depends on the reduced density ma- 
trix at previous times t' . In the stationary limit, how- 
ever, all matrix elements of p re d become time indepen- 
dent (Pg(t) = PJg{t') = P* 1 ), and with the additional 

definition W Xl x } = f £ dt' W XlX Ut,t'), Eq. © sim- 

X2X2 J ~°° X2X 2 V ' ^ 

plifies to 

= i p £ = -\ - *w p ti + E Kli p i -(8) 

xi,x 2 

The stationary charge current through lead r is given by 

^= e t wi xi p i>' 0) 

where the matrix elements of the current transition rates 
W Ir are directly obtained by multiplying the correspond- 
ing matrix elements of W with the net transported charge 
from lead r to the island. As we consider the stationary 
state, the charge on the island is conserved. Hence we 
are able to define the current flowing through the device 
as / = II = — Ir- For a detailed derivation of this di- 
agrammatic language and its rules for the calculation of 
diagrams, we refer to Refs. E^^S 

For a metallic single-electron transistor with a dense 
level spectrum on the island, the dimension of the re- 
duced density matrix and, thus, the number of kinetic 



equations determining them is huge. A drastic simplifi- 
cation of the problem is achieved by getting rid of the 
off-diagonal matrix elements of the reduced density ma- 
trix on the right-hand side of the kinetic equations. Since 
the tunneling Hamiltonian conserves charge and is diag- 
onal in the transverse channel index, only those matrix 
elements P x , of states |x) = \{ni a u}} and \x') = \{n f lav }) 
with the same number of island electrons in any channel 
v need to be considered, ni au = n\ av . A further 
simplification relies on the assumption that there is a fast, 
spin-independent energy relaxation within the island, i.e., 
the time scale of energy-relaxation r er is smaller than the 
dwell time ra w (which, in turn, is smaller than the intrin- 
sic spin-flip time r s f in order to sustain spin imbalance 
on the island). This is reasonable due to the fact that 
in contrast to quantum dots, metallic islands accommo- 
date a large number of electrons and hence exhibit many 
relaxation channels. As a consequence, any coherent su- 
perposition between states with different occupations of 
the orbital levels / are destroyed and only coherences be- 
tween different spin states is kept. This leads to the more 
restrictive condition S ^ J(7 ni av = n' l(7V . Since the spec- 
trum on the island is spin degenerate, energy differences 
E x — E x f, that would appear on the right-hand side of 
the kinetic equation, vanish. 

To remove the remaining coherent superpositions of 
different spin states we choose the spin quantization axis 
on the island along the direction of the accumulated 
spin and, furthermore, assume that the steady-state spin 
structure of the island is rotationally invariant about the 
quantization axis. The latter assumption neglects any 
anisotropics of quadrupole and higher moments within 
the plane perpendicular to the dipole moment. Eventu- 
ally, we conclude ni au = n' l(7V for the proper choice of 
the spin quantization axis, i.e., only diagonal matrix el- 
ements P x = P* enter the right-hand side of the kinetic 
equations in the stationary limit, 

° = i p £ = Y, w m p x> (io) 

x 

and similar for the current. 

The assumed fast, spin-independent energy relaxation 
within the island does not only destroy coherences of 
states with different occupations of the orbital levels 
it also leads to a thermal equilibrium among all states 
that are connected by this relaxation, i.e., states with 
given numbers and N± of spin ^ and 4- electrons 
on the island. The individual occupation for level I 
with spin a and channel index v under the condition 
that the island contains and N± electrons with spin 
t and 4, respectively, can be expressed by the condi- 
tional probabilities F(lcriy\N^^ N^). Beenakker first in- 
troduced these probability functions to discuss resonant 
tunneling through a quantum dot coupled to two electron 
reservoirs^ while Barnas et al. used them in the con- 
text of collinear single-electron spin- valve transistors^ 
In thermal equilibrium, F{lau\N^, N±) reduces to the 



Fermi function with a spin-dependent chemical poten- 
tial, f(e l - fJi a (N a )) with f(E) = l/[exp(^) + 1], The 
spin-dependent chemical potential fi a (N a ) is determined 
by the condition 



N a = Y,fh-^(Na)} 



(11) 



The kernels in Eq. ([T0|) depend on the initial 

state x- But most of the information about the individual 
occupations of the levels contained in x are irrelevant for 
the evaluation of the kernel matrix element. What mat- 
ters is the occupation and the energy of only those island 
levels that are involved in the tunneling processes that 
take place in the transition described by the considered 
kernel. In this paper, we restrict ourselves to first order in 
the tunnel-coupling strength. Each kernel is, then, a sum 
over contributions for which only one level I and trans- 
verse channel index v is involved (we remind that coher- 
ent superpositions of different levels I or transverse chan- 
nels v do not appear). In the notation \x) = Kxz^;^)), 
which makes the occupation of level I and channel index 
v explicit, the non- vanishing kernel matrix elements can 
be written in the form i.e., the part xiv 

is not changed during the transition. The value of the 
kernel depends on the total island charge N = + N± 
via the charging-energy contribution to the Hamiltonian 
and the occupation and energy of only the selected level 
I and transverse channel index v. This occupation, how- 
ever, is fully determined by the spin-dependent chemical 
potential [i a (N a ), which, in turn, depends on and N±. 

For an island with a dense level spectrum, the depen- 
dence of the chemical potential fi a (N a ) on N a is rather 
weak: the change for each added electron with spin a is 
of the order of the mean level spacing Ae and, thus, much 
smaller than the energy scales ksT^eV relevant for trans- 
port. Therefore, we use fi a (N a ) « fi a independent of N a 
in the following. This has the consequence that the value 
of the kernel matrix element depends only on the 

total island charge N (of state x-> i- e -? N = N° with 

= T,iv(x\4av c i™\x)), the energy and occupation of 
the island level I with channel index v that is involved 
in the transition. This dependence can be cast in the 
notation W^$(h v -> N), where the occupation of the lev- 
els and channels different from I or v does not appear 
explicitly. 

The independent degrees of freedom governed by the 
kinetic equations are the probabilities to find TV electrons 
on the island, 
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(12) 



as well as the three components of the total island spin, 
or, equivalent ly, the magnitude and the direction of the 
island spin. The latter enter the right-hand side of the 
kinetic equations through the spin splitting A/i = /^ — /i^ 
of the chemical potential, S = ftp/A/i/2, and the angles 
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FIG. 3: (Color online) Energy scheme of single-electron spin- 
valve transistor. The island is described as two independent 
electron reservoirs related to the spin species a G {t>i}- The 
chemical potentials /jll,r describe the applied bias voltage and 
A/x = n^—ni determines the spin accumulation on the central 
electrode. The charge states N are tunable via gate voltage 
and €f denotes the Fermi energy. 



a and /? appearing in the tunnel matrix elements. We 
remark that, for a constant density of states on the island, 
the average chemical potential \i = (/^ + /xjJ/2 does not 
depend on the amplitude of the accumulated spin S. For 
each independent degree of freedom we need one kinetic 
equation. These are provided by the kinetic equations 
for the following operators A: the projector \N) (N\ = 
J2 X \x)(x\^n,n x on the state with total island charge N 

and the total island spin S = (ft/2) Y^ioo'v c \<ju°° °' c i°'vi 
where a is the vector of Pauli spin matrices. The kinetic 
equations are given by 



o = -p N = J2s N , Nx w^; Px 



At 

A , A 



(13) 



°=d^ )= £ (X2\S\xi)W^,P x ,. (14) 



Introducing all the simplifications on the right-hand side, 
as discussed above, and making use of the notation 
f+(ei) = f(ei - fjL a ) and f~(ei) = 1 - f(e t - fi a ) leads to 

0= o^ = E<o°G^^)/f(Q)A"(Q)^ 



W^(l, u,N- l)/ t -(eO/ 4 -(Q)^-i 
Y,WrA^N)f+{e l )f^{e l )P N 

lucr 

W^(l, v,N — l)/+(c,)/f (cO^V-i 
WZiil, v,N+ l)f+(e l )f+(e l )P N+1 
£ W^(l, u, NV+Wf+MPN , (15) 
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for the probabilities Pn- For the z component of the 
island spin we obtain 



( xw; (;,^iv)/ t -(e / )4-(e i )^ 

r 

hJ2aW™(l, v, N)f+{ti)fr{ti)P N 

Nlva 

KJ2°W^(1, v, N)f+(ei)f+ (ei)P N ,(16) 



Nlua 



Nlua 



where a contributes to sums with a + sign for a =t and 
a — sign for a =|. The kinetic equations for the x— and 
^/-components of the island spin, expressed by the raising 
and lowering operators (S^ = (S x ± iS y ), are 



= 



■.-(s+)=hJ2wf°(i,v,N)f^(e l )f^(e l )p N 

Nlv 
Nlva 



N 



Nlv 

and the kinetic equation for (S~) is the same but replac- 
ing Wf$(l,v,N) with Wl$(l,v,N) on the right-hand 
side. 

The explicit values of W$*$(l,v,N) are given in Ap- 
pendix [A] We plug them in and replace the summation 
over the island levels I by an integration over an energy 
cj, multiplied by the density of states of the given spin pj . 
In the present paper, we assume the central electrode to 
be nonmagnetic, hence the island density of states is spin 
independent, pj" = pi/ 2. The occurring integrals can be 
carried out and we finally get the following master equa- 
tions: 



d p 



ir^2 [a+ CT (Ajvr_i)Pjvr_i + a ra (A N )P N+1 

ra 

-a+ a (A N )P N - cv^Ajv-OPjv] , (18) 



— (S z ) = irh * K*( a n) - ^(Ajv-i)] Pn- 



Nr 



Here, we used the definition of the island rate functions 



at(E) := ±a° ro 



E - (fl r - p a ] 


) 


exp 


_ k B T 


- 1 



(19) 



where a® 



PiNc 
2-kH 



T 1 ! is the dimensionless conductance 



of lead r for island spin <j, N c the number of trans- 
verse channels, and E is the energy of the tunneling elec- 
tron. The equations (fl8]) are linearly dependent. Hence, 
to solve the master equations we need an additional 
equation that is provided by the normalization condi- 
tion ^ N Pn = 1. For collinear single-electron spin- valve 
transistors ((/> G {0,7r}), the spin accumulation naturally 
points in a direction parallel to the lead magnetization 



directions n r , this directly defines the angles a and f3. 
Therefore, the equations in Eq. (fl8|) are sufficient to de- 
scribe such a system but not the noncollinear case. 

In the derivation of Eqs. ([T6]) - (fT8]) . we have used a 
specific coordinate system for the spin. In the chosen 
coordinate system, the x and y components of S vanish. 
The final version of the kinetic equations for the spin 
can, however, be written in a coordinate-free represen- 
tation. First, we proceed in a manner analogous to the 
paragraph above, i.e., the kernels are plugged in and the 
summation over I is replaced by an integration over en- 
ergy uj. Starting from the master equations of (S ± ), see 
Eq. (fT7|h and by using the kinetic equation of (S z ), see 
Eq. JTHJ), one finally obtains the following coordinate- free 
representation of the master equation for the island spin 
S: 



d(S) /d(S) 



dt 



dt 



d(S) 
dt 



rcl 



d(S) 
dt 



• (20) 



Accumulation, relaxation, and rotation processes deter- 
mine the time evolution of S. The contribution, which 
builds up an average spin, reads 



Ms>\ 

V dt J, 



p r -— [n r + (n r • n s )n s \ 



Nra 



[a- (Ajv-i) - a+ CT (A N )] P N . (21) 



Transitions to a charge state N by tunneling of electrons 
from the leads onto the island and vice versa are de- 
scribed by this term and lead to a polarization of the 
island. On the one hand, tunneling processes of both 
leads r accumulate a spin in the direction of the respec- 
tive polarization n r . But on the other hand, due to 
the macroscopic spin of the central electrode, there is 
an additional accumulation contribution in the direction 
of S. In equilibrium (V = 0), there is no current flow- 
ing through the central electrode and the equation above 
becomes (^jf ) = 0. Hence, without any applied bias 
voltage, there is no spin accumulated on the island, i.e. 
S = 0. In contrast to the accumulation contribution, the 
relaxation term 

V al /rel Nrv J 

x[/-(w + Ajv-i) - + Ajy)] ,(22) 

causes a decay of the island spin, which is characterized 
by s(cj). The vector 

sM = y[/tM-AH]fis (23) 

corresponds to the energy-dependent spin density in the 
central electrode. The influence of the exchange field 
between ferromagnetic lead r and island is represented 
by the third contribution of Eq. ([20]) : 



V dt ) 



h 



£ |d W sHxBL 



:H- (24) 
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Here, we used the dimensionless magnetic moment of 
electrons g, the Bohr magneton and the definition 
of the exchange field between island and lead r: 



f 



N 

frW) 



V/V-l 



where the prime at the integral denotes a principal value 
integral. The exchange field can be interpreted as a 
many-body interaction effect. In the limit of noninter- 
acting island electrons (Ec = 0) it vanishes independent 
of gate voltage. The contributions B£ xc (cj) act on the 
spin like an applied external magnetic field, which points 
in the polarization direction of lead r. This results in a 
precession of the accumulated spin out of the (ulAr) 
plane. Equation ([25]) is similar to the expression of an 
exchange field existent between a single-level quantum 
dot and ferromagnetic leads, which was first introduced 
by Braun et al. in the context of a quantum-dot spin 
valve4£ But in contrast to the quantum-dot spin valve, 
there is a net spin accumulated on the island in all possi- 
ble charge states P/v- Hence the total exchange field is a 
composition of all charge-state contributions. Addition- 
ally, the integration over uj in Eq. ([24]) results from the 
continuous density of island states. 

Having solved the master equations for P/v and (S), the 
stationary charge current through lead r may be calcu- 
lated. We start from I r = ^ %%/ {W Ir )** f P x > . Execut- 
ing the same procedure as in the previous considerations 
of the master equations one gets 



I r = [ a ta( A N) ~ a r<T (A N -i)] Pi 

No 



N 



(26) 



IV. RESULTS 

We consider electric transport both in the linear and 
the nonlinear regimes. For this, we calculate the accu- 
mulated spin S and the (linear and nonlinear) electric 
current as a function of the system parameters such as 
gate voltage Vg, bias voltage V, angle <p between source 
and drain's magnetization direction, and the degree of 
polarization p. An explicit focus in the discussion is put 
on the influence of the exchange field on the transport 
characteristics. 



A. Linear-response regime 

We start by considering the single-electron spin-valve 
transistor in the linear-response regime (eV <C fceT, Ec). 
All quantities are periodic in the gate voltage Vq with pe- 
riodicity e/Cc- For low temperatures, at most two charge 
states Nq and Nq + 1 have a non-vanishing occupation 
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FIG. 4: (Color online) Magnitude of exchange field in linear 
response as a function of applied gate voltage for different 
angles between the lead polarization directions. Parameters 
Pl — Pr — 0.3, Yl — Yr — r/2, and Ec — lO/c^T were 
chosen. 



probability. All analytic formulas for the linear-response 
regime are derived for this limit. With increasing tem- 
perature, more charge states may become occupied. The 
results plotted in the figures are always calculated taking 
into account all relevant charge states. 

To describe the linear-response regime, the system of 
master equations Eqs. (fl8|) and ([20]) is expanded up to 
first order in the transport voltage V. The term associ- 
ated with rotation of the accumulated island spin simpli- 
fies to 



d(S) 
dt 



lin 



(27) 



where we used the expression of the exchange field in 
linear response 



B 
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(28) 



The magnitude B l ™ c = |B^ C | is plotted for different an- 
gles as a function of the gate voltage in Fig. HI The 
exchange field vanishes at the symmetry point between 
two resonances (middle of Coulomb valley), i.e., at in- 
teger values of CgVgI^^ in accordance to the case of 
the quantum-dot spin valve. At half-integer values of 
CgVg/g, where the two charge states Nq and Nq + 1 are 
degenerate (An = 0), the exchange field vanishes due 
to a cancellation of the Nq and Nq + 1 contributions of 
Bg XC (w). This is in contrast to the case of a quantum- 
dot spin valve, where the exchange field is maximal at 
the resonance position for transport. To understand this 
difference we consider the resonance between charge state 
(empty dot) and 1 (singly occupied dot) of the single- 
level quantum dot. Here, the symmetry is broken due 
to the fact that on the one hand the state 1 can be vir- 
tually excited to the doubly-occupied state 2 but on the 
other hand there is no corresponding virtual excitation 
of the state in a lower charge state. Hence the different 
contributions do not cancel each other. 



7 



n L xn R 



FIG. 5: (Color online) Scheme of accumulated spin S in linear 
response regime. The spin precesses within a plane defined 
by fix, — fiR and nz, x fiR. The spin precession is accompanied 
by a decreased magnitude of accumulated spin. 



In the following, we assume symmetric polarizations 
(Pl = PR = p) and tunnel-coupling strengths (Yl = 
Tr = r/2). In this case, the total exchange field points 
in the + n^-direction and rotates the spin out of the 
(ul^ur) plane, see Fig. [5] The island spin S acquires 
a component perpendicular to the (ul^r) plane. We 
obtain a lm = and 



/3 lin 0) = -arctan 



V T7V C A 



No 



k B T 



.(29) 



In Fig. [6(a), the rotation angle is plotted as a function 
of the gate voltage. Sign changes appear both at inte- 
ger and at half-integer values of CcVc/e, which reflects 
the sign change of the exchange field at these positions. 
Away from these points, the angle f3 hn is finite. Close 
to but not exactly at integer values of CcVc/e and at 
low temperatures (Ec ~ lO/c^T), the angle approaches 
±7r/2. There, a small variation Vg induces a strong re- 
orientation of the spin S, as illustrated by the insets of 
Fig. [6(a). This high sensitivity can be quantified by the 
slope of f3 hn with respect to Vg- For low temperature 
(k B T <C Ec), we find that the slope is proportional to 
k B Tp cos | exp (Ec /k B T)/ Eg, which increases exponen- 
tially with decreasing temperature. As a consequence, 
the gate- voltage range close to integer values of CcVc/e 
is ideal for manipulating the direction of the island spin. 

The absolute value 5, which is determined by the accu- 
mulation and relaxation term of the spin master equation 
reads 



> lin W = ^^sin^cos/3 lin ((/>). 



(30) 



It is proportional to the cosine of /3 lm , hence the magni- 
tude decreases with increasing angle f3 hn and vanishes for 
pim _ see Fig. [5] Thus, the effect of the exchange 

field is not just an rotation of S but also a reduction of 
its magnitude S. In Fig. [6(b), # lin (7r/2) is plotted for 
various temperatures. We see that for Ec < 3k B T the 
structure is smeared out and the absolute value of spin 
is nearly constant due to a suppression of the Coulomb 
blockade. 

After considering the accumulated spin, we now an- 
alyze the linear conductance G hn = (dI/dV)\v=o of 



CO 



I 
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FIG. 6: (Color online) (a) Angle f3 hn and (b) spin accu- 
mulation S hn for different temperatures as a function of 
gate voltage. The chosen parameters are pl = Pr = 0.5, 
F L = F R = T/2, and tt/2. 



the single-electron spin- valve transistor. The conduc- 
tance, normalized to the conductance of a single-electron 
spin-valve transistor with parallel lead magnetizations 
G lin (0) = e 2 N c p I TA N J(8k B Tsmh^), can be ex- 
pressed in terms of the rotation angle /3 lm (</>) as 



G lm (0) 
G lin (0) 



= 1 



p 2 sin 2 § 



1 + tan 2 /3 lin (0) 



(31) 



In Fig. [7(a), the conductance G hn ((j) = tt/2) is plotted 
as a function of gate voltage. Since it is an even function 
of charging energy Aat , the conductance is symmetric 
with respect to the resonance points, given by half-integer 
values of CcVc/e. Away from the resonances, the con- 
ductance is suppressed due to Coulomb blockade. The 
difference between the solid and dashed lines illustrates 
the influence of the exchange field. The dashed lines are 
obtained by manually setting B l ™ c to zero. The exchange 
field rotates the accumulated spin out of its blocking po- 
sition. Hence, it increases the conductance of the single- 
electron spin-valve transistor except for the symmetry 
point at resonance. As a consequence, for high polariza- 
tions, p > 0.9, the maximum at resonance can even turn 
into a local minimum. 

Consideration of the conductance as a function of the 
magnetization angle see Fig. [7(b), shows that the 
exchange field reduces the spin- valve effect. We empha- 
size that, though the solid and dashed lines coincide for 
4> G {0, 7r, 27r}, the exchange field does not vanish in these 
points. In fact, it points parallel to the accumulated spin 
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FIG. 7: (Color online) Linear conductance in units of e 2 TpiN c 
with (solid) and without (dashed) the effect of the exchange 
field: (a) Vb-dependence in case of <j) — tt/2 and (b) <j>- 
dependence with applied gate voltage CgVg — 4e/10. For 
both plots the parameters pl = Pr = 0.9, Yl = Fr = T/2, 
and Ec — 15k bT were chosen. 




FIG. 8: (Color online) (a) Current I through single-electron 
spin- valve transistor in units of eTpiN c for different polar- 
izations p, (b) the corresponding island spin accumulation S, 
and (c) the relevant occupation probabilities Pn in the case 
p — 0.3. For all three plots, the parameters pl — Pr — P, 
Y L = F R = T/2, = tt/2, C g Vg = 3e/10, and E c = 50k B T 
were chosen. 



and, therefore, Bg£ c does not rotate S out of the (n^n^)- 
plane. 



B. Nonlinear- response regime 

We now turn to the nonlinear-response regime, eV > 
ksT^Ec- The transport voltage V is symmetrically 
applied to both leads, i.e., the chemical potentials are 
fiL = eV/2 and {ir = —eV/2. The current through 
the single-electron spin- valve transistor as a function of 
bias voltage for different polarizations p is shown in Fig. 
Ufa). Here, the angle between the lead magnetizations 
is chosen to be <p = tt/2. Due to the larger spin ac- 
cumulation for higher lead polarizations, the current re- 
duces with increasing p. At low bias voltages, the trans- 
port through the system is blocked until eV exceeds the 
Coulomb blockade threshold 2Atv and both spin reser- 
voirs of the next charging state (No + 1) enter the trans- 
port window. We introduce the notation En ](J for the 
excitation energy of the a channel of charge state N. 
From this definition it follows that Ejvo+i^ = 2Aat . 
By further increasing of the source-drain voltage, more 
and more levels of the continuous spectrum contribute 
to transport, and the current increases continuously (in 
contrast to a stepwise increase for quantum dot with a 
discrete level spectrum). We note that we assumed here 
symmetric tunnel couplings, Tl = Tr = T/2. By an 
asymmetric choice of the tunnel couplings, the Coulomb 
step at eV = 2An and also the steps of the other charg- 



ing states that enter for higher V can be made much 
more pronounced i 82 ^ 83 In Fig. [He), the corresponding 
occupation probabilities of the relevant charging states 
(N e {N - 2, N - 1, iV , N + 1}) are plotted as a func- 
tion of the bias voltage (for p = 0.3). The probability 
to find the island in a given charge state N decreases 
for states with higher excitation energies. In general, the 
voltage that is necessary to excite the island state (TV, a) 
is determined by the equation eV = \2An + aAfi(V)\. 
Here, A/x(V) is the island level splitting for the respective 
bias voltage and a contributes with a factor +1 for the 
t reservoir and -1 for the 4- reservoir. It is proportional 
to the spin accumulation 5(0), which is plotted in Fig. [8] 
(b). By means of A/x(V) the transition voltages can be 
determined self consistently. We emphasize that due to a 
finite island spin accumulation, either the spin "f- or spin 
| reservoir enters the transport window at the transition 
voltages of the occupation probabilities, see Fig. 03 

An analysis of the bias- voltage dependence of the is- 
land spin [see Fig. [8] (b)] yields that S can be decom- 
posed into two components, an oscillating one and a 
component that monotonously increases ] 37 ! 38 Since the 
former is suppressed for higher p, let us consider the 
graph of p = 0.1 to explain this behavior. Our theory 
describes sequential tunneling, hence in the Coulomb- 
blockade regime, spin accumulation is exponentially sup- 
pressed. For eV > 2Aat , both spin reservoirs of the 
charging state No + 1 contribute to transport and S in- 
creases with increasing V until the t reservoir of charging 
state No — I enters the transport window. Now, there is 
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an additional channel for the t electrons to leave the is- 
land and S decreases until the corresponding 4- reservoir 
is reached. Subsequently, the spin accumulation increases 
until 4 electrons can tunnel into charging state No + 2. 
A reduction of S follows, then the state t, No + 2 can be 
occupied and the spin accumulation increases again. By 
further increasing of V this processes of depletion and 
filling repeat. To understand the reduced influence of 
the oscillating component for higher p, we first assume 
antiparallely polarized leads (j) = u. In this case, the tun- 
neling rates of the processes, which lead to a decrease of 
S, are proportional to 1 — p, hence they occur less for 
higher p. Since there is a more complex but similar be- 
havior of the tunneling rates for finite noncollinear angles 
(/>, we can conclude that, in general, the oscillations are 
weaker for higher polarizations. 

The finite spin accumulation on the island results in a 
nonzero TMR. Due to the noncollinear system we define 
an angular-dependent magnetoresistance: 



TMR(0) 



1(0) -W) 
1(0) 



(32) 



with I ((f)) being the current through the single-electron 
spin-valve transistor in which the magnetization direc- 
tions of the leads enclose the angle (j). In Fig. [9j the 
bias- and gate-voltage dependence of the TMR is plot- 
ted for different polarizations and an angle (j) = tt/2. 
In both cases, the TMR shows an oscillatory behavior 
originating from the bias/gate dependence of the spin 
accumulation^ 32 1 36 Naturally, the magnetoresistance is 
more pronounced for higher p. For lower polarizations, 
we observe sign changes of TMR caused by the oscillating 
spin accumulation, 37 We want to emphasize that for ob- 
taining reliable results in the Coulomb-blockade regime 
[marked by the grey area in Fig. [9fa)], one has to take co- 
tunneling processes into account, which we neglect in our 
theory. Next, we focus on the analysis of the influence of 
the exchange field on the TMR, which is illustrated by the 
dashed-dotted lines in Fig. [9f b). To obtain the lines rep- 
resenting the case of an absent exchange field, we manu- 
ally set Bg XC = in Eq. ([20]) . For gate voltages represent- 
ing a vanishing exchange field (CqVgI^ £ {0,1/2,1}), 
naturally, the two graphs coincide. Comparison of the 
TMR gate dependence for different polarizations shows, 
that the exchange field has a stronger effect for higher 
p. By affecting the accumulated spin the exchange field 
decreases the TMR. 

In order to identify the threshold voltages at which 
new transport channels open, it is convenient to study 
the second derivative of the current d 2 I/dV 2 , see Fig. 
[TOl The positions of the peaks directly represent the al- 
ready discussed excitation energies of the relevant charg- 
ing states. All charging states induce two peaks, one 
for each spin channel (t?i)- An exception is the first 
peak (No + 1), which does not split into two due to 
the exponentially suppressed spin accumulation at the 
relevant voltage, see Fig. [5Jb). The vertical lines in 
the figure mark the voltages at which the charge states 
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FIG. 9: (Color online) TMR of a noncollinear setup (cj) — 
7r/2) for different polarizations p. Plot (a) shows the bias 
dependence for CgVg = 3e/10 and (b) the gate dependence 
for eV = 2Ec- The dashed-dotted lines represent the cases 
where the exchange field is manually set to Bg XC = 0. For 
both plots, the remaining parameters were chosen to be pl = 
PR =P,T L =r R = T/2, and E c = 50k B T. 
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10: (Color online) Second derivative of the current 
d 2 I/dV 2 in units of e 3 TpiN c for two different polarizations 
p and <j) — tt/2. The excitation energies of the charge/spin 
states represented by the peaks strongly depend on the lead 
polarization p. Parameters pl — PR = p, Tz, = Tr = T/2, 
CgVg = 3e/10, and Ec = 50k bT were chosen. 



(TV e {TVo - 1, N + 1, N + 2}) would enter the transport 
window if one neglects spin accumulation. Due to the 
fact that the accumulated island spin strongly depends 
on the lead polarizations, the excitation energies Ej^-a- 
are also sensitive to a variation of see Fig. [TTJ The 
described behavior suggests that in experiments, a mea- 
surement of d 2 I/dV 2 as a function of bias voltage can 
be used as a convenient tool to determine the spin split- 
ting of the chemical potential on the island and, thus, 
the degree of polarization of the leads. 
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FIG. 11: (Color online) Polarization dependence of the excita- 
tion energies En-,o- normalized to Ec- The chosen parameters 
are p L = p R = p, = tt/2, Y l = F R = T/2, C G V G = 3e/10, 
and E c = 50k B T. 



V. CONCLUSION 

We have investigated electronic transport through a 
noncollinear single-electron spin- valve transistor realized 
by a metallic island weakly coupled to two ferromag- 
netic leads. Employing a diagrammatic real-time tech- 
nique, we performed a systematic perturbation expan- 
sion in the tunnel-coupling strength. We found that an 
interaction-induced exchange field appears as a conse- 
quence of a spin-dependent tunnel coupling to the leads. 
In the linear-response regime, this exchange field weakens 
the spin- valve effect and increases the linear conductance. 
It, furthermore, leads to a high sensitivity of orientation 
of the accumulated island spin on variation of the gate 
voltage Vg, which is convenient for controlled manipu- 
lation of the island spin. Finally, we suggest the second 
derivative of the current d 2 I/dV 2 in the nonlinear regime 
as a tool to determine the spin splitting of the chemical 
potential on the island, which, in turn, depends on the de- 
gree of spin polarization in the leads of the single-electron 
spin- valve transistor. 
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(for the kinetic equations of P/v and S z ) or off-diagonal 
(for the kinetic equations of 5' ± ). Explicit calculation 
yields 



W^{l,v,N) = 


fl 

rs 


h A N ), 


(Al) 


W§Z(l,u,N) = 


rs 


h Ajv_i), 


(A2) 


W$Z(l,v,N) = 


rs 


h Ajv), 


(A3) 


WZ$(l,u,N) = 


rs 


1- Ajv_i), 


(A4) 



where f^~(E) = f r (E) describes tunneling of electrons 
from lead r into the island and fr(E) = 1 — f r (E) tun- 
neling out of the island into lead r. Hence the elements 
above describe processes that change the occupation of 
the central electrode. These parts of the kernel are re- 
lated to the diagonal matrix elements where the total 
island charge remains constant via 



Wgg(l,v,N) =-J2w^(l,u,N), 



(A5) 



WZZ{l,u,N) =-W§Z(l,v,N) - Wti(l,v,N),{A6) 
Wjj{Uv,N) =-Y,WZ d d {l,v,N). (A7) 



The matrix elements with off-diagonal final states can 
also be divided into those where the total island charge 
changes during the transition, 

W?°(l,v,N) = — J2p r sV s r :V^f+(e l + A N ), (A8) 



2tt 



Y,P r s V sa V s-Jr^l + AjV-l),(A9) 



and those where the total island charge remains constant, 

w^{l,v,N) = -^^y;>^/(iv), (Aio) 

rs 

w;t{i,v,N) = ^x>:v£v£j>), (ah) 



Appendix A: Kernel elements of master equations 

In the calculation of the matrix elements of the kernel 
W to first order in the tunnel coupling only one island 
level I and channel number v is involved. Its value de- 
pends on the energy and the occupation of this level as 
well as on the total island charge. As discussed in the 
main text, we need matrix elements that describe tran- 
sitions from a state that is diagonal in the occupation of 
the considered island level to a state that may be diagonal 



with the integral expression 



ei + A N -i -L0 + M+ 



ei + An — cu — z0+ 



(A12) 



which represents transition processes of island-charge 
states, while Re I(N) exclusively describes virtual charge 
transfer. 
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